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ABSTRACT 

We have compared numerical simulations to observations for the nearby (< 40 Mpc) 
groups of galaxies (jHuchra &: GelleHll98 2' and'R ameha et al.ll2002f) . The group identi- 
ficatio n is carried out using a group- finding algorithm developed bv lHuchra fc Gelleil 
(|l982f ). Using cosmological N-body simulation code with the ACDM cosmology, we 
show that the dynamical properties of groups of galaxies identified from the simulation 
data are, in general, in a moderate, within 2a, agreement with the observational cat- 
alogues of groups of galaxies. As simulations offer more dynamical information than 
observations, we used the N-body simulation data to calculate whether the nearby 
groups of galaxies are gravitationally bound objects by using their virial ratio. We 
show that in a ACDM cosmology about 20 per cent of nearby groups of galaxies, 
identified by the same algorithm as in the case of observations, are not bound, but 
merely groups in a visual sense. This is quite significant, specifically because estima- 
tions of group masses in observations are often based on an assumption that groups of 
galaxies found by the friends-of-friends algorithm are gravitationally bound objects. 
Simulations with different resolutions show the same results. We also show how the 
fraction of gravitationally unbound groups varies when the apparent magnitude limit 
of the sample and the value of the cosmological constant A is changed. In general, a 
larger value of the Q\ generates slightly more unbound groups. 

Key words: methods: numerical - galaxies: clusters - galaxies: haloes - cosmology: 
dark matter - large-scale structure of Universe. 
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1 INTRODUCTION 

Small groups of galaxies are the most common galaxy 
associations and contai n ^ 50 p e r cen t of all galaxies 
in the universe (e.g. iHolmberd Il950l: iHumason et al.l 



19561: iHuchra fc Gellej Il982l : iGeller fc Huchral Il983l : 
Nolthenius fc White! Il987h . The study of galaxy groups is 
a very interesting area of research because these density 
fluctuations lie between galaxies and clusters of galaxies 
and may provide important clues to galaxy formation. 
Small groups of galaxies are also important cosmological 
indicators of the distribution, and properties, of dark 
matter in the universe. 

The dynamics of the nearby groups of galaxies and the 
Local Group has provided a unique challenge to cosmological 
models in the p ast. The quiescence of the local peculiar ve- 
locity field (e.g. Ide Vaucouleurj[l958l : ISandage fc TammannI 
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1 19751 : ISandaed Il999l and lEkholm et al] I2001I ) was a long- 
standing puzzle that presented a challenge for the mod- 
els of structure formation. The velocity field within 5h~^ 
Mpc of the Local Group is extremely 'cold', the dis- 
persi on is only ~ 50 — 60 km s~^ (|Teerikorpi et al.] 
I2OO5I and references therein). The ACDM cosmology and 
dark energy has solved this p r oblem and it has been 
shown by e.g.lKlvpin et all (120031), iMaccio et all l|2005l ') and 
IPeirani fc de Freitas Pachecol ( ioOo ) with constrained sim- 
ulations that the ACDM cosmology can produce the small 
values of the velocity dispersion. Today, the question of the 
virialization of groups of galaxies and the fraction of grav- 
itationally bound systems provides a new challenge for the 
cosmological models and grouping algorithms. 

Over the last two decades, cosmological simulations 
have proven to be an invaluable tool in testing theoreti- 
cal models in the non-linear regime. The standard approach 
is to assume a cosmological model and to use the appro- 
priate power spectrum of the primordial perturbations to 
construct a random realisation of the density field within a 
given simulation volume. The evolution of the initial den- 
sity field is then followed by using an N-body simulation 
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code, and the results in the simulation box (viewed from 
outside) are compared with observational data. A compar- 
ison of the simulations with observational data is typically 
done in a statistical manner. The statistical approach works 
well if there is a statistically representative sample of objects 
with well-understood selection effects for both the observed 
universe and the simulations. 

Given a set of observed galaxies with their positions 
in the sky and their redshifts, the task of a group-finder 
is to return sets of galaxies that most likely represent 
true gravitationally bound structures. Some contamina- 
tion is always expected due to selection effects in obser- 
vations. This study uses one of the most popular group- 
finders: the friends-of-friends (FOF) algorithm. FOF has 
been used wi dely for identifying groups of galaxies in the red- 
shift surveys (Huchra & Geller 1982';'Geller & Huchra"l98d; 
Nolthcnius & White 1987; RamcUa ct al 1989 ; Moore ct aL, 
"l993| : iRameUa et al.lll997i:lRainella et al]|2002l : iTucker et"al] 
20od and iGiuricin et aLmOOof T and is now a standard ap- 



proach. 

The identification of the group members has, in gen- 
eral, been based on a su bjective selection of dat a. In order 
to remove this difficultv. iHuchra fc Gelleil l|l982l ) (hereafter 
HG82) developed a method of identifying groups of galaxies 
from the observations. It has been usually thought that the 
FOF based algorithms would produce groupings which are 
mostly gravitationally bound when the number of galaxies i n 
a group exceeds five members (see e.g. Ramella et al.ll2003 ). 
A few studies (e.g. ICarlberg et al.l 1200 ll ) have been made 
where the free parameters of the FOF algorithm have been 
optimised to avoid spurious groups with interlopers. Studies 
of the FOF algorithms have concluded that the choice of the 
free parameters depends on the case that is studied, and no 
singular best choice of the parameters can be made. Unfor- 
tunately, none of the observational methods do truly answer 
the question whether groups of galaxies are gravitationally 
bound objects. There are some estimates for the fraction of 
how ma ny groups found by th e FOF algorithm are spurious 
(see e.g. iRamella et"aLll2002l ). but these estimates are not 
based on the physical properties of the groups. 

Cosmological N-body simulations include all the nec- 
essary information for finding out whether a given object 
is gravitationally bound or not. lAceves fc Velazquez! (|2002l ) 
studied small galaxy groups with N-body simulations and 
made conclus i ons a bout the virialization of the groups. 
iDiaferio et al.l (119941 ) studied compact groups of galaxies 
with N-body simulations and made remarks about the frac- 
tion of bound groups and chance alignment systems. Groups 
of galaxies, and their dynamical properties generated by the 
FOF algorithm, have been studied earlier with N-body sim- 



l y ! 

ulations (see e.g. Nolthenius et al.ll 19971: Diaferio etahll 19991 : 



iMerchan fc Zandivarej|2002l : ICasagrande fc Diaferio! l2006l ). 
Most of these earlier studies have taken advantage of the 
constrained simulations and have used models which are 
different from the currently popular ACDM cosmology. In 
these constrained, or mock catalogue, simulations the sim- 
ulated space has been compared directly with the observa- 
tions of redshift space and the groups have been categorized 
as spurious if the algorithm has failed to find similar groups 
as in real space. Instead, we have studied the virial ratio 
of the groups of galaxies produced by the FOF algorithm 
first described in HG82. We will show that ~ 20 per cent of 



the groups of galaxies found by the FOF algorithm are not 
real gravitationally bound groups, but spurious. Our results 
agree roughly with the pre vious results a n d es timates but do 
not confirm the claim bv I Ramella et af] (!l989„) that groups 
with more than four members are gravitationally bound. 

For an 'observer' placed in a specific location, selecting 
a similar environment between observational data and cos- 
mological simulations might be problematic. The simplest 
way is to choose an 'observation' point within a simulation 
box by certain criteria. It is arg ued fe.g. !Klvpin et al.!!2003! 
and references therein) that it is not clear what 'similar en- 
vironment' actually means and that simply placing the ob- 
server at some specific point would resolve the issue. In this 
paper, we show that this is a useful approach, as we are 
not comparing simulations to observations directly but with 
statistics. 

The main purpose of this work is to study if groups of 
galaxies found by the HG82 algorithm are bound, and how 
the fraction of bound groups depend on the chosen magni- 
tude limit and cosmological model. We will show our findings 
with different apparent magnitude limits and for different 
cosmological models. We also show that there is no signif- 
icant correlation between the crossing time of a group and 
its virial ratio. 

This paper is organized as follows. In section[2]we review 
the method used in the identification of the group members 
in the observations. A brief discussion of the differences be- 
tween observations and simulations is given in section (2] In 
section O we discuss briefly the virial ratio, used for deter- 
mining whether a group of galaxies is gravitationally bound. 
Section |4] discusses the simulations we use for the analysis. 
In sections O and |6l we present our results and discuss the 
findings. Discussion of the probability functions of gravita- 
tionally unbound groups is done in section [71 Finally, we 
summarize our results in section |8] 



2 A REVIEW OF THE GROUP-FINDING 
ALGORITHM 

In observations there are generally three basic pieces of in- 
formation available for the study of the galaxy distribution: 
the position, the magnitude and the redshift of each galaxy. 
Although the magnitude is important as a measure of the ob- 
ject's visibility, it is usually a poor criterion for group mem- 
bership. The method used for creating a group catalogue in 
HG82 can be summed up in two criteria: the projected sep- 
aration and the velocity difference. The FOF algorithm is 
described in greater detail by the original authors in HG82. 

The grouping method begins with a selection of an ob- 
ject, which has not been previously assigned to any of the 
existing groups. After choosing the object the next step is 
to search for companions with the projected separation D\2 
smaller or equal to the separation Dl- 



Di2 = 2sin(^-^ — ^ DL(Vi,l^2,mi,m2) , 
where the mean cosmological expansion velocity is 



(1) 



(2) 



and the velocity difference V12 is smaller or equal to the 
velocity Vl'- 
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\Vi-V2\^VL(Vi,V2,mi,m2) , 



(3) 



where Vi and V2 refer to the velocities (redshifts) of the 
galaxy and its companion, mi and m2 are their magnitudes, 
and 6 is their angular separation in the sky. If no companions 
are found, the galaxy is entered in a list of isolated galaxies. 
All companions found are added to the list of group mem- 
bers. The surroundings of each companion are then searched 
by using the same method used in the first place to find com- 
panions. This process is repeated until no further members 
are found. 

There are a variety of prescriptions for Dl and Vl ■ We 
adopt the method used in HG82, and assume that the lumi- 
nosity function is independent of distance and position and 
that at larger distances only the fainter galaxies are missing. 
For each pair we take: 



Dl = Do 



(4) 



where the integration limits can be calculated from equa- 
tions: 



Miim = miim - 25 - 5 log(DF) 
and 



M12 — ran 



25-51og(\/) 



(5) 



(6) 



and where <E>(M) is the differential galaxy luminosity func- 
tion for the sample, and Dq is the projected separation in 
Mpc chosen at some fiducial distance Df- In this paper we 
adopt constants Do = 0.63 Mpc and Dp — 10 Mpc to be 
same as in HG82. The effect of varying Do has been stud- 
ied e.g. by iRamella et all l|l989l ) where the effects are also 
explained. 

The limiting velocity difference is scaled in the same 
way as the distance D^: 



Vl^Vo 



J*2>™ HM)dM ^ 



(7) 



where the fiducial value is Vo ~ 400 k m s ^ and the integra - 
tion limits are as above (Eq. [5] and |6} . iRamella et"aLl IjlQS^ ) 
varied Vo and concluded that the results are not sensitive to 
the choice of Vo- This is probably related to the geometry 
of the large-scale structure. Frederic (1995a,b) argues that 
the optimal choice of Do and Vo depends on the purpose for 
which groups are being identified. Similar claims has been 
made in pape rs where the FO F algorithm has been opti- 
mized (see e.g. I Eke et al 1 l2004l and lBerlind et al.ll2006l ). Be- 
cause of this, we also show some results when Do = 0.37 Mpc 
and Vo = 200 km s~^ are adopted. 

Note that the scaling law in Equations |4] and [7] 
has been questioned by many authors. Specifically, re- 
pla cing the power -1/3 by -1/ 2 (see the argument in 
e.g. iNolthenius fc White! HqStI and iGourgoulhon et allll992l ) 
drastically reduces the correlation between the redshift and 
the velocity dispersion observed in the HG82 group cata- 
logue. However, part of this correlation is related to a selec- 
tion effect rather than to the grouping algorithm, because 
groups with low velocity dispersion usually have few bright 
galaxies and so they can be seen only at low redshift. In this 



paper, we use the equations mentioned above for consistency 

with the HG82 catalogue. 

For simplicity, we use the ISchechteil (|l974 ) luminosity 
function: 



$(M) = -$*lnlO 
5 



exp 



-10!^ 



(8) 



where M is the absolute magnitude of the object. We adopt 
the parameter values of a = —1.02, M* — —19.06, and 
$* = 0.0277 comparable to IIG82. For comparison, we 
use the galaxy luminosity function values of a = —1.15, 
M* = —19.84, and $* = 0.0172, wh ich were derived fro m 
the Millennium Galaxy Catalogue bv [Driver et al.l l|2007l ). 

Our chosen values of constants and parameters needed 
for the FOF algorithm are exactly the same as in HG82 
for consistency. We do show selected results with more re- 
cent values of constants and parameters if these results differ 
from the results produced with the IIG82 values. Through- 
out this paper we adopt the parameterized Hubble constant 
Ho — lOOh km s~^ Mpc"^ with h = 1.0 for comparison 
with HG82 when a definite value of the Hubble constant is 
needed. We do not consider dust extinction in the analysis 
of the simulations in this paper. 



3 A REVIEW OF THE VIRIAL THEOREM 

In the simplest case when we have a two body system with 
total mass M = mi + m2 and V is the relative speed of 
the components, the kinetic energy T of the system (in the 
centre-of-mass coordinate system) and its gravitational po- 
tential energy U (taken positive) are: 



^ 1 mim2 2 
2 M 



U = 



Gm\m2 
R 



(9) 
(10) 



where R is the size of the system and G is the gravitational 
constant. These equations are related by a simple relation: 



U = 2T 



(11) 



However, the above relation holds only for isolated self- 
gravitating systems when the system is in equilibrium. 

In general, groups of galaxies (and dark matter haloe^ 
contain more than two members. Therefore a general- 
ized method for calculating the kinetic and the poten- 
tial energies is needed. W e adopt a method described by 
[Ch ernm fc Mikkolal l|l99lD . In general, the kinetic energy 
may be written: 



and the potential energy: 



(12) 



(13) 



^ We will use the term 'halo' from now on to refer to virialized 
clumps of dark matter in the simulation and reserve 'galaxy' for 
the real observational data. 
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where rrii and rrij are the masses of the two galaxies, V i 
and V j are their velocities, and Ri,-j is the distance between 
them. 

We use these general equations to get the total kinetic 
energy of a group of haloes and compare it to its total po- 
tential energy. If the group of haloes does not fulfill the cri- 
terion: 

T-f/<0, (14) 

it is entered into a list of unbound groups. The above crite- 
rion is equal to the virial ratio: 

J < 1.0 , (15) 

which we use throughout this paper as a criterion for dis- 
criminating between bound and unbound groups. 



Table 1. Details of the simulations analyzed in this paper. 





L 




Zi 




^ res 


Nh 


0.73 


40 


256^ 


47.9600 


2.86 X 10** 


1.8 


6700 


0.73 


80 


256^ 


38.7099 


2.29 X 10^ 


7.3 


9301 


0.90 


80 


256^ 


31.7858 


8.50 X 10* 


7.3 


4937 


0.00 


80 


256^ 


72.8767 


8.47 X lO'' 


3.7 


7919 



Note: Q\ specifies the value of the cosmological constant, L is the 
size of the simulation box in one dimension in h~^Mpc, Np is the 
number of dark matter particles, Zi is the initial redshift, rrires 
is the mass resolution in h~^MQ, Fres is the force resolution 
in h~^kpc, and N/^ is the total number of dark matter haloes 
identified from the simulation. 



4 DESCRIPTION OF THE COSMOLOGICAL 
SIMULATIONS 

4.1 Background 

We present results from four simulations, performed by 
the cosmological N-body simulation code AMIGA (Adap- 
tive Mesh Investigations of Galaxy Assembly). The former 
version of AMIGA was known as MLAPM (for details see 
iKnebe et ai]|200lj ). For the first two runs we adopt the cur- 
rently popular flat low-density cosmological model ACDM 
with h = 1.0, D.dm = 0.27, Qa = 0.73 and as = 0.83, with 
two different resolutions. Both simulations were made with 
256^ dark matter particles. The high resolution simulation 
began at the initial redshift of Zi = 47.96 while the low- 
resolution simulation was initiated at redshift Zi — 38.71. 
The volume employed in the high resolution simulation 
was {40h~^ Mpc)^ and {80h~^ Mpc)^ in the low resolu- 
tion simulation corresponding to the mass resolutions of 
2.86 X 10*/i"^ Mq and 2.29 x lO^h'^ Mq, respectively 
The force resolution for the high resolution simulation is 
1.8h~^ kpc and for the low-resolution 7.3h~^ kpc. 

For the third and the fourth simulations we adopt dif- 
ferent cosmological models. These simulations were also per- 
formed with 256^ dark matter particles but with different 
values of the cosmological constant A. The total density of 
the universe was kept equal to the critical density (SI = 1.0). 
For the third simulation we adopt h = 1.0, Qdm ~ 0.1, 
JIa ~ 0.9 and as = 0.83. For the fourth simulation we adopt 
h = 1.0, Qdm = 1.0, D.A = 0.0 and erg = 0.84 (see Table [ij. 

The high resolution SIa = 0.73 simulation was used to 
understand the effects of limited resolution in N-body simu- 
lations. During this work, we found some differences between 
results of the high and low-resolution Qa =0.73 simulations. 
These differences are clearly visible when the group abun- 
dances are studied (see Figures 1 — 4). 

4.2 Halo finder and identification of tiie dark 
matter haloes 

Our simulations only follow the evolution of the dark mat- 
ter particles via gravitational interaction. It is expected 
that baryons condense and form galaxies at the centres of 
dark matter haloes. AMIGA N-body simulation code comes 
with a halo finding algorithm called MHF (MLAPM's Halo 



Finder. [Gill et al]|2004h . For our purpose of analyzing the 
nearby groups of haloes we used MHF. 

The general goal of a halo finder, such as MHF, is to 
identify gravitationally bound objects. MHF essentially uses 
the adaptive grids of the AMIGA to locate the haloes and 
the satellites of the host haloes, namely subhaloes. The ad- 
vantage of reconstructing the grids to locate haloes is that 
they follow the density field with the exact accuracy of the 
simulation code and therefore no scaling length i s required. 
For m ore detailed description of the MHF see I Gill et aU 
(! 2004l ). 

The minimum number of particles in a halo was set 
to 10. This corresponds to a halo mass ~ 3 x 10^h~^ Mq 
and ~ 2 x 10^°h~^ Mq for the high and the low-resolution 
simulation, respectively. A low value of the minimum num- 
ber of particles in a halo ensures that even with a limited 
resolution, large and massive haloes are split into lighter 
subhaloes. For large and massive (~ 10^^/i~^ Mq) haloes, 
subhaloes represent visual galaxies, masses ^ 10^^h~^ Mq. 
In a typical case, the total massfraction in subhaloes is 
~ 10 per cent, and only these are visible in our 'mock' cat- 
alogue. In our 'mock' catalogue, the median of individual 
dark matter haloes mass is ~ 1.6 X lO^^h'^ Mq when the 
Qa = 0.73 model and the low-resolution is adopted. The 
first and the third quartiles are: ~ 7.2 x 10^^ Mq and 
~ 5.0 x lO^^fe"^ Mq, respectively. 

The AMIGA and its halo finder calculate automatically 
certain properties (e.g. position, mass, velocity etc.) of the 
dark matter haloes. These properties were used when the 
FOF algorithm was applied to generate the catalogues of 
groups of dark matter haloes. Subhaloes were included in 
our data as our purpose is to study if the groups of galaxies 
(dark matter haloes) are gravitationally bound objects. The 
results did not change substantially when the subhaloes of 
the more massive haloes were excluded from the analysis. 
This result is due to the small number of subhaloes in our 
low-resolution simulation. The results of the high resolution 
simulation show no significant difference if subhaloes were 
excluded because subhaloes have a relatively small mass. 
Due to their small masses, subhaloes are not visible at the 
observation point, when the apparent magnitude limit of 
13.2 is adopted. 
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5 STATISTICAL PROPERTIES OF THE 
NEARBY GROUPS OF GALAXIES: 
COMPARING SIMULATIONS WITH 
OBSERVATIONS 

5.1 Selection of the nearby groups of haloes 

A Total of ten catalogues were generated, corresponding to 
10 different 'observers', for each simulation, with different 
apparent magnitude limits. Ten observers are used to pro- 
duce enough groups to give good statistics. Note, however, 
that since all ten catalogues are constructed from the same 
parent simulation, the scatter between statistics estimated 
from them might underestimate the true sampling variance. 

All observation points were chosen with the following 
criteria: 

(i) observation point is > 15ft~^ Mpc from the edge of 
the simulation box 

(ii) a massive (~ 5 x lQ^^h~^ Mq) Virgo -type halo is 
located within a distance of ~ 20h~^ Mpc. 

We did not restrict the local < 10h~^ Mpc environment of 
the observation points by any criteria. Although it is not 
clear if choosing an observation point simply by the two for- 
mer criteria resolves the environment issue, we believe this 
to be strict enough for the statistical study of the virial 
ratio of groups. These criteria are justified for the statis- 
tical study, as we did not perceive a significant difference 
between the observation points in the low resolution sim- 
ulation. Small differences between observation points were 
observed when the high-resolution simulation was studied, 
as it only includes two Virgo -type haloes. When the low- 
resolution simulation was studied, the location of a massive 
(Virgo -type) halo did not have any significant effect to our 
results. 

The simulation data do not directly give the luminosity 
or the absolute magnitude of the dark matter haloes, which 
are needed when we are mimicking observational conditions. 
We use halos' virial mass M^ir to obtain its luminosity. To 
obtain the luminosity of an object in the blue band we use 
the relation proposed by I Vale fc Ostrikerl (|2004| ): 



L{Myir 



■ 5.7 X W^h-^LQ 



q + M,l 



(p-r) 



1/s 



where Mn is defined: 



Mil 



10"/i-iM, 







(16) 



(17) 



For the free parameters of the mass- luminosity function, val- 
ues of P = 4. 0, q — 0.57, r — 0.28, a nd s = 0.23 were adopted 
('OEuri"200^). It has been shown bv'Coor av fc Milosavlievia 
(2005) that the relation between the mass of a dark matter 
halo and its luminosity is not as straightforward as presented 
above. For our purposes, as the luminosity of a dark matter 
halo is used only to determine whether a halo is visible from 
the observation point, the above relation should be satisfac- 
tory. For this work we do not adopt more complex methods 
such as actual distributions for the mass-luminosity relation. 

After the luminosity L of the halo is known we obtain 
the apparent magnitude of the halo in the blue band from 
the equation: 



ma 



-2.51ogi 



+ 51ogi 



d 



1 Mpc 



-h25 ,(18) 



where d is the distance from the observation po int and th e 
magnitude of the sun in blue band Mq^ = 5.47 (|Coxll2000l l. 
As seen from Equation [TS] we do not include dust extinction 
in our study, as its effect in a statistical study like ours would 
be negligible. 

The method described above allows us to use the ap- 
parent magnitude limit mum ~ 13.2 as adopted in HG82. 
Group catalogues of this study are also generated with dif- 
ferent magnitude limits in order to understand the effects 
of the magnitude limit in magnitude-limited samples. Un- 
less explicitly noted, all haloes and groups referred to are 
from the simulatio ns; real gr oups of galaxies from HG82 and 
UZC-SSRS2 (|Ramella et al.|[200a) are denoted as such. 



5.2 Comparison parameters 

We begin by calculating the velocity dispersion of a 
group. In general, the velocity dispersion of a group is de- 
fined as: 



\ 



Nh 



N„ 
i = l 



(19) 



where Nh is the number of haloes (or galaxies) in a group, Vi 
is the radial velocity of the ith halo (or galaxy) and < > 
is the mean group radial velocity. 

The second comparison parameter is the mean pairwise 
separation Rp which is a measure of the size of a group. It 
can be defined as: 



„ S <VR> . 

Rp = — — sm 

TTiin 



1 



N„{Nh 



j<i i—l 



(20) 



where < vr > is the mean group radial velocity. Ho is the 
Hubble constant, and 9ij is the angular separation of the ith 
and jth group members. Other two comparison parameters 
are the total group mass and the virial crossing time (in 
units of the Hubble time Hn^) which can be defined as: 



3Rh 
53/2cr, 



(21) 



where ct„ is the velocity dispersion and Rh is the mean har- 
monic radius: 



„ TV <VR > . II 

^« = — H^^^M 2 



Nh{N„ - 1) j 



= 1 j>i 



where all the variables are defined as above. 

In observations, the group masses can be estimated in 
various methods. In the HG82 and the UZC-SSRS2 cata- 
logues the total mass of a group is estimated with a simple 
relation: 



Mobs = 6.96 X 10 a„i?j/M0 



(23) 



where is the velocity dispersion, and Rh is the mean 
harmonic radius of the group as defined above. We use this 
simple relation to determine the 'observable' mass of a group 
when we are comparing the masses of the simulated groups 
to the real observed groups such as HG82 and UZC-SSRS2. 



6 Sami-Matias Niemi et al. 



However, when the total mass of a group is used to scale 
the properties of groups (in Section [6)l, we calculate the to- 
tal mass of a group of dark matter haloes as a sum of the 
member halo masses, namely the 'true' mass of the group. 
In case of the subhaloes, the diffuse dark matter is included 
in the main halo's mass. In general, we do not include the 
diffuse dark matter within a given distance from the group 
centre, as it would be troublesome to choose an appropriate 
distance. We do not consider this error to be meaningful, as 
the diffuse dark matter does not substantially give rise to 
the mean density of a simulation. 



5.3 Comparison with observations 

Comparison between simulations and observations are done 
using a Kolmogorov-Smirnov (K-S) test. The null hypothesis 
Hnuii of the K-S test is that the two distributions are alike 
and are drawn from the same population distribution func- 
tion. Results of the K-S tests are presented as significance 
levels (value of the Q function) for the null hypothesis. Cor- 
relation between two variables is proved or disproved with 
the use of the linear correlation coefficient r. In general, 
the significance level of O.OOf is adopted when the correla- 
tion between two variables is determined. For correlations 
we also present a probability P{r) of observing a value of 
the correlation coefficient greater than r for a sample of 
observations with N — 2 degrees of freedom. 

We begin the comparison of our simulations to obser- 
vations by using the parameters presented in the previous 
subsection. A direct comparison between HG82 and simula- 
tions is possible as the magnitude limit (13.2) and the depth 
of the catalogues {cz < 4, 000 km s~^) are comparable. Com- 
parisons with the more recent group catalogue UZC-SSRS2 
is also done. The UZC-SSRS2 catalogue has the magnitude 
limit of 15.5 and only galaxies with cz < 15, 000 km s"'^ has 
been considered. These differences make the direct compari- 
son of the UZC-SSRS2 with the simulations less conclusive. 
We also compare our ten observation points with each other 
but do not find significant difference between them. This 
justifies our method of choosing the observational points as 
stated before. 

In Figures [1] — [5] the abundance of groups is scaled to 
the volume of a sample as the distributions depend strongly 
on selection and volume effects. However, as there is no 'to- 
tal' volume of a galaxy sample in magnitude-limited group 
catalogues, we wei ght each group accordin g to its distance 
l|Moore et al.ll 19931 and lDiaferio et al.ll 19991 ). As we only con- 
sider groups with three or more members, we can identify a 
group only when its third-brightest galaxy has an absolute 
magnitude 



Mi < mii 



25 



.< cz > 



(24) 



where < cz > is the mean velocity of the group. Mi deter- 
mines the radius czi of a sphere within which we could have 
identified this group. 

We calculate the comoving volume sampled by a group 
with the following equation: 



3 \Ho 



3zi / Q 
2 V 2 



of the group, c is the speed of light, and f2 is the cosmological 
density parameter, taken as 1.0. Each group of galaxies (or 
dark matter haloes) contributes with a weight of ^P"^ to the 
total abundance of groups. We include all galaxies with cz > 
500 km s~^. This lower cut-off avoids including faint objects 
that are close to the observation point as these groups could 
contain galaxies fainter than real magnitude-limited surveys. 
Therefore we consider only groups with < cz > larger than 
500 km s"^ in mock, HG82, and UZC-SSRS2 catalogues. 

Figures [1] — 14] show that the ACDM simulations are, in 
general, in a moderate agreement with observations when 
the resolution effects of the N-body simulations are taken 
into account. Our ACDM simulations are within 2a from the 
UZC-SSRS2 catalogue and within Scr from the HG82 cata- 
logue. Figures [1] — 14] are all from the SIa = 0.73 simulations 
when the apparent magnitude limit has been set to 13.2, 
comparable to HG82. The error bars in Figures [l] — |4] are the 
standard deviation between ten observation points. Unless 
explicitly noted, the simulations referred to are = 0.73 
simulations, other models are denoted as such. 



5.3.1 Velocity dispersion 

Figure [1] shows that the cosmological Qa = 0.73 
model can pro duce velocity disper s ions similar to obse rva- 
tions (see also Klvpin e t al. 200^4 iMaccio et al .' 2005 and 
IPeirani fc de Freitas Pacheco 20061). Fig ure [T] agrees with 



(25) 



where to is the solid angle of the catalogue, Zi is the redshift 



results by Casagrande fc Diaferid (l2006li (their Figure 14) 
even though Casagrande fc Diaferid ! 20061 ) considered only 
groups with > 5 members. Our low-resolution simulation 
produces roughly the right number density of groups when 
the high (> 100 km s~^) velocity dispersions are considered 
and the comparison is carried out against more recent ob- 
servations (UZC-SSRS2). However, due to the limited mass 
resolution, the low-resolution simulation lacks a significant 
number of groups when the abundance of groups with ve- 
locity dispersions < 100 km s~^ is studied. Because of this 
discrepancy the applied K-S test fails: Q ~ 10~^ (against 
the HG82) and Q ~ 10"^ (against the UZC-SSRS2). Even- 
though the K-S test fails, the low-resolution simulation is 
within 3(7 from the HG82. When the high-resolution simu- 
lation is considered we get roughly the same number density 
of groups as in observations. However, the high resolution 
simulation lacks groups with velocity dispersions > 500 km 
s~^. This can be explained by the small volume of the high- 
resolution simulation. Even with this discrepancy, the ap- 
plied K-S tests are approved with the significance levels of 
0.02. When observations (HG82 and UZC-SSRS2) are com- 
pared against each other, the applied K-S test is approved 
at level of 0.34. For detailed significance levels of the K-S 
tests see Table |31 



5.3.2 Mass 

Figure [2] shows that the cosmological ACDM model can pro- 
duce 'observable' masses (Eq. I23|l similar to observations 
when the resolution effects of the simulations are consid- 
ered. The low-resolution simulation can produce the same 
number density of groups as in the UZC-SSRS2 catalogue 
when massive (log(Group Mass M"^) > 13.5) groups are 
considered. Both the UZC-SSRS2 catalogue and the low- 
resolution mock catalogue has an excess of massive groups 
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Figure 1. The cumulative number density of velocity dispersion 
Ov for galaxy groups. Simulation data are from the = 0.73 
simulations, and it is averaged over the ensemble of ten mock 
catalogues. The error bars are Itr errors and are only shown for 
the low-resolution Oa = 0.73 simulation for clarity. The error 
bars for other data have similar size. 

if comparison is carried out against the HG82 catalogue. 
When less massive groups (log(Group Mass M'^) < 13.0) 
are studied, the low-resolution simulation has a number den- 
sity of groups which is over Scr lower than in the UZC-SSRS2 
catalogue. The resolution effect is clearly visible in Figure 
[2] when the high resolution simulation is studied, as it can 
produce about the right number density of groups when less 
massive (log(Group Mass M"^) < 13.0) groups are consid- 
ered. The high-resolution simulation is less than Icr away 
from the HG82 and within 2cr from the UZC-SSRS2 even 
when groups with log(Group Mass Mg^) < 11.0 are consid- 
ered. The applied K-S test is approved {Q ~ 0.58) only when 
the high-resolution simulation is compared to the HG82 cat- 
alogue. In all other cases the K-S test fails. For numerical 
details of the K-S test, see Table |3] 

As the 'observable' mass of a group depends strongly on 
the groups' velocity dispersion (see Equation I23|l we made 
another comparison between group abundances by mass. If 
we use the 'true' mass of a group instead of an 'observ- 
able' mass, other differences arise. There is no substantial 
difference between the plots of 'observable' and 'true' mass 
when comparing the abundances of lighter (log (Group Mass 
M-i) < 14.0) groups. However, our simulations do not pro- 
duce a single group with a 'true' mass > 5 x 10^''/i~^Mq. 
The low and high-resolution SI a = 0.73 simulations show the 
same cut-off, thus the lack of massive groups is not a res- 
olution effect. However, the small volume of our simulation 
boxes explains the lack of massive groups in simulations and 
the existence of groups with large 'observed' mass is due to 
projection effects, which are not realiably taken into account 
in Eq. ^ 

5.3.3 Size 

The cosmological Q.a ~ 0.73 model can produce groups of 
haloes which are similar in size to observed groups. Note, 
however, that we do not compare our simulations to the 
UZC-SSRS2 catalogue, as it doesn't contain the informa- 
tion about the pairwise separtion of groups. Our high- 
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logio(Group Mass / IVIass of the Sun) 



Figure 2. Group abundance by 'observable' mass of the groups. 
Simulation data are from the Oa = 0.73 simulations, and it is 
averaged over the ensemble of ten mock catalogues. The error bars 
are Icr errors and are only shown for the low-resolution Ca = 0.73 
simulation for clarity. The error bars for other data have similar 
size. 

resolution simulation produces about the right number den- 
sity of groups when small (log(Rp) < —0.4) groups are con- 
sidered and the error is well within la. The high- resolution 
simulation seems to produce an excess of groups when in- 
termediate size (log(Rp) £ [—0.3,0.4]) groups are consid- 
ered. However, as the only comparison observation is the 
HG82 catalogue this excess might not be as large as in Fig- 
ure |3] as other comparisons (Figures [1] and [2]) show that the 
HG82 and the UZC-SSRS2 catalogues differ quite signifi- 
cantly from each other. 

When the low-resolution simulation is studied, we ob- 
serve this same excess when larger (log(Rp) > 0.1) groups 
are considered. Because of these discrepancies the applied 
K-S test fails in both cases, with Q ^ 10~^ and Q ~ 10~^ 
for the high and the low-resolution simulations, respectively. 
Even though the K-S test fails, the simulated mock cat- 
alogues of group abundances by mean pairwise separation 
are mostly within 2a. 

5.3.4 Crossing time 

Figure U shows that the cosmological S7a ~ 0.73 model can 
produce groups with crossing times similar to the HG82 
observations. The low-resolution simulation produces the 
number density of groups with small crossing times, which 
is a lot lower than observed. However, this discrepancy is 
due to limited resolution, as the high-resolution simula- 
tion produces a lot more groups with small crossing times. 
The high-resolution simulation produces roughly the right 
number density of groups when the crossing time of the 
group is studied. Some differences are observed when larger 
(log(tc) > —0.6) crossing times are studied. Both simula- 
tions produce a higher number density than observed. For 
low-resolution simulation this excess is not large as the num- 
ber density is within 2a. Because of the discrepancies visible 
in Figure |4l the applied K-S test fails in both cases. For nu- 
merical details, see Table |3] When more recent values of the 
Schechter luminosity function are adopted, somewhat lower 
crossing times are observed in general. However, more re- 
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Table 2. The correlation between crossing time and virial ratio 
when different values of Oa and the apparent magnitude limits 
are adopted. 



0.2 
Logio( Rp)(h"'Mpc) 

Figure 3. The cumulative distribution of mean pairwise sep- 
aration Rp for galaxy groups. Simulation data are from the 
= 0.73 simulations, and it is averaged over the ensemble of ten 
mock catalogues. The error bars are Itr errors and are only shown 
for the low-resolution Ha = 0.73 simulation for clarity. The error 
bars for other data have similar size. 
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Figure 4. The cumulative distribution of crossing time tc (in unit 
of the Hubble time) for galaxy groups. Simulation data are from 
the r^A = 0.73 simulations, and it is averaged over the ensemble 
of ten mock catalogues. The error bars are la errors and are only 
shown for the low-resolution Ha = 0.73 simulation for clarity. 
The error bars for other data have similar size. 



cent values of the Schechter luminosity function do not give 
a better agreement, and the K-S test fails. The significance 
levels of the K-S tests are ~ 10~* and ~ 10~^, respectively. 

Our simulations with the FOF algorithm do not (with 
a few exceptions) contain groups with crossing time larger 
than one Hubble time. For the high-resolution simulation, 
the median value of the crossing time is 0.14,11^^. The me- 
dian value of the crossing time for the HG82 catalogue 
and for the low-resolution simulation is Q.l^H^^ . Small 
crossing times suggests that groups of galaxies have had 
time to virialize a nd these grou ps shoul d be gravitation- 
ally bound (see e.g-ICott fc Turner 1977: Tucker et al-llioool : 



lAceves fc Velazauej|2002l and lPlionis et al.ll2006l ). We stud- 



f^A 




r 


a 


P{r) 


0.00 


13.2 


-0.009 


> 0.05 


0.33 


0.00 


20.0 


0.005 


> 0.05 


0.36 


0.73^ 


13.2 


-0.016 


> 0.05 


0.26 


0.73^ 


13.2 


-0.016 


> 0.05 


0.29 


0.73^ 


20.0 


-0.016 


> 0.05 


0.09 


0.90 


13.2 


-0.012 


> 0.05 


0.42 


0.90 


20.0 


-0.022 


> 0.05 


0.09 



ied the correlation between crossing time and virial ratio 
and did not find any significant relation between these two 
variables. The linear correlation coefficient of —0.01 suggests 



Note: Oa specifies the value of the cosmological constant (^ = 
high-resolution and ^ = low-resolution simulation), is the 

apparent magnitude limit of the search, r is the value of the lin- 
ear correlation coefficient, a is the significance level, and P{r) is 
the probability of observing a value of the correlation coefficient 
greater than r. 

that there is no correlation between crossing time and virial 
ratio, in the low-resolution simulation when the apparent 
magnitude limit of 13.2 is adopted. This correlation is not 
significant at level of 0.05 and P{r) ~ 0.29. There is no 
significant correlation between these two variables when dif- 
ferent values of Qa, resolutions or the apparent magnitude 
limits are adopted (see Table [2]). The lack of correlation 
be tween virial ratio a nd crossing time (see similar results 
m iDiaferio etmiiggal ) calls into question the crossing time 
as an estimator of gravitationally bound systems which is 
widely accepted in observations. 



5. 3. 5 Richness 

The Q,A ~ 0.73 model can produce groups comparable to 
observations when the number of members in a group is 
studied. The abundance of 'rich' groups is roughly the same 
in the low-resolution simulation as in the observations. How- 
ever, the low-resolution simulation cannot produce as many 
'poor' groups as is observed. The lack of 'poor' (< 4) and 
the excess of 'intermediate' (g [6,40]) groups is the reason 
why the K-S test fails {Q ~ 10~^). The agreement is even 
worse {Q ~ 10"^) when more recent values of Schechter lu- 
minosity function are adopted. These values (a = —1.15, 
M* = —19.84, and $* = 0.0172) produce a large number of 
'poor' groups and a lack of 'rich' groups. The difference be- 
tween the low-resolution simulation and observations is due 
to the limited resolution in the simulations. When the high- 
resolution simulation is used, the agreement to observations, 
especially to UZC-SSRS2, is better (Q ~ 10"^). 



5. 3. 6 Influence of dark energy 

There are big differences between different cosmological 
models when dynamical properties of groups of dark mat- 
ter haloes are studied. Figure [5] shows the impact of dark 
energy on the formation of galaxy groups. It is clear that 
the f^A = 0.0 simulation over produces groups with high ve- 
locity dispersions. The excess is over 3a if the comparison 
is carried out to the HG82 catalogue. A smaller discrep- 
ancy is observed when the comparison is carried out to the 
UZC-SSRS2 catalogue. The small number density of small 
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A = 0.00 
1a 

A = 0.90 
1a ■ 
HGS2 
UZC-SSRS2 




Table 3. Comparison of HG82, UZC-SSRS2, and simulations 
when the apparent magnitude limit of 13.2 is adopted. 



1.4 1.6 1.8 2 2.2 2.4 2.6 2.8 3 3.2 
Logio( ) (km s"') 

Figure 5. The cumulative distribution of velocity dispersion cr^ 
for galaxy groups from two different simulations with two differ- 
ent values of the cosmological constant f2A when the apparent 
magnitude limit of 13.2 has been adopted. The simulation data 
are averaged over the ensemble of ten mock catalogues. For com- 
parison also the observational data from the HG82 and the UZC- 
SSRS2 catalogue have been plotted. The standard = 0.73 
simulations are shown in Figure [T] 



f2A 




HG82 


UZC-SSRS2 


HG82 vs. UZC-SSRS2 


0.00 


0"u 


10-5 


lo-'' 


0.34 


0.73^ 


(Tv 


0.02 


0.02 


0.34 


0.73^ 


o"„ 


lo-'^ 


lo-'' 


0.34 


0.90 


cr„ 


10-^ 


10-7 


0.01 


0.00 


Ma 


10-5 


10-s 


10-3 


0.73^ 


Mg 


0.58 


10-5 


10-3 


0.73^ 


Mg 


lo-** 


10-^ 


10-3 


0.90 


Mg 


10-*^ 


10-7 


10-3 


0.73^ 


Rp 


10-3 






0.73^ 


Rp 


10-5 






0.73^ 


tc 


10-" 






0.73^ 


tc 


10-5 







Note: Significance levels of the K-S test for the null hypothesis 
that observations and the simulations = high-resolution and 
^ = low resolution) are alike and are drawn from the same parent 
population (HG82 and UZC-SSRS2 columns). Significance levels 
of the K-S test for the null hypothesis that the HG82 and the 
UZC-SSRS2 group catalogue are alike and are drawn from the 
same parent population {HG82 vs. UZC-SSRS2 column). 



velocity dispersion groups can be explained by the resolu- 
tion effect which is also visible in Figure [1] Because of these 
discrepancies the K-S tests fails. When the Q.a = 0.90 sim- 
ulation is studied, a qualitatively better agreement is ob- 
served, especially when the comparison is carried out to the 
HG82 catalogue. The great difference in the number density 
of small velocity dispersion groups can be explained by the 
resolution effect and the fact that the 0,^ = 0.90 simula- 
tion has relatively small number of groups. Because of the 
great discrepancy in the number density of small velocity 
dispersion groups, the applied K-S test fails in both cases. 

The Q.A =0.0 cosmology produces more massive groups 
with greater velocity dispersion than the ACDM cosmology. 
Larger number of massive groups can partially be explained 
by the somewhat lower mass resolution in the f^A = 0.0 sim- 
ulation. However, the excess of massive groups is most likely 
due to Equation 1231 which we use to obtain the 'observable' 
mass of a group, which depends strongly on the velocity dis- 
persion of the group. For numerical details of K-S tests, see 
Table 13 



5. 3. 7 Median values and other properties 

The median values of the group properties are presented in 
Table m In general, our simulations seem to produce groups 
which median value of the velocity dispersion and the group 
mass is greater than in observations. In simulations, groups 
have also a greater median value for the mean pairwise sep- 
aration than in the IIG82 sample. The Oa = 0.73 simula- 
tions have the median value of velocity dispersions which 
are close to observations, even though they are somewhat 
higher. In general, median values of the group properties are 
in a moderate agreement with the res ults of similar studie s 



Table 5. Fractions of isolated galaxies, binary galaxies, and 
groups. 



(e.g. Casagrande fc Diafer"iol 20061 and lDiaferio et ai]|l999l ) 



Source 


'^lim 


Groups (%) 


Binaries (%) 


Isolated {%) 


Simu 


13.2 


38.1 ± 0.5 


19.7 ±0.1 


42.2 ±0.5 


Simu 


20.0 


41.2 ± 0.2 


17.6 ± 0.1 


41.2 ±0.2 


HG82 


13.2 


60.0 


14.0 


26.0 


LEDA 


14.0 


33.3 ±0.5 


17.1 ± 0.3 


49.6 ±0.4 



ICasagrande fc Diaferid l|2006l ) found larger values for the 



median velocity dispersions and group masses, but they con- 



Note: Source refers to the sample (Simu = the low-resolution 
f^A = 0.73 simulation), is the apparent magnitude limit 

(in m^jQ) except for LEDA in By), Groups (%) is the fraction 
of galaxies in the groups. Binaries (%) is the fraction of galax- 
ies forming double systems, and Isolated (%) is the fraction of 
galaxies which are not classified into any group or double system. 
Poisson error limits have been calculated for the samples. 

sidered only groups with > 5 members, which most likely 
makes the median values of groups somewhat higher. 

The fractions of isolated galaxies, binary galaxies and 
groups of galaxies were also studied. If we compare our 
results with HG82, a significant difference in the fraction 
of iso lated galaxies is noticed. Comparison to the LEDA 
(|Giud icc 1999) ca talogue shows a be tter fit. For more de- 
tails, see Table [S] iTucker et al.l (|2000D listed a large fraction 
of galaxies in groups from different group catalogues. These 
results and comparisons are not shown in this paper due 
to the different magnitude limits and grouping algorithms 
adopted in those observations. We may state that, in general, 
simulations show similar results to observations (excluding 
HG82), with regard to the fractions of groups, binaries and 
isolated galaxies. 

We also made an attempt to study discord ant redshifts 
in compact groups observed e.g. by iSulentici (|1984) ) and 
iGirardi et al.l (Il992[). This effect has been studied by several 
autho rs ( see e.g. "iBvrd fc Vahonenlll985l . IVahonen fc BvrdI 
Il986l andE ovmo fc Hicksonlll997l ) who have come up with 
different explanations. According to these authors, apparent 
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Table 4. Weighted quartiles of the group properties. 





= 0.00 


Ha = 0.73^ 


Ua = 0.73^ 


CIa = 0.90 


HG82 


UZC-SSRS2 


log Mg 
Rp 

tc 


178/307/403 
13.6/14.0/14.4 
0.85/1.38/1.57 
0.09/0.14/0.25 


135/180/295 
12.7/13.2/13.5 
0.81/1.21/1.47 
0.07/0.14/0.19 


103/160/298 
12.8/13.6/14.0 
0.80/0.99/1.24 
0.10/0.19/0.27 


124/209/271 
12.7/13.1/13.5 
0.72/1.01/1.20 
0.09/0.17/0.23 


60/135/155 
12.0/12.4/13.2 
0.44/0.54/1.00 
0.10/0.19/0.27 


60/130/178 
12.1/12.8/13.3 



Note: Quantities a^, log Afg, Rp, and tc are in units of km s h ^Mq, h ^ Mpc, and Hubbies, respectively. 



discordant redshifts arise when groups are not virialized and 
their central galaxies are incorrectly identified. Our findings 
are not conclusive as we did not have any exact method to 
identify which dark matter haloes might represent observ- 
able spiral galaxies. We did not observe any significant asym- 
metry in the radial velocities of the groups and neither this 
asymmetry was seen in the groups, which were misidentified 
(so that the brightest member is not the dominant member). 
No significant difference for the radial velocity asymmetry 
was discovered between bound and unbound groups. 



6 GRAVITATIONALLY BOUND GROUPS 

Gravitationally bound groups are determined by using the 
criterion (virial ratio, Eq. I15|l presented in Section [S] This 
method of computing the gravitational potential well of 
a group does assume that the group is isolated. This is 
not strictly true as each group is embedded in the large- 
scale matter distribution, which might have an effect to the 
threshold 1.0 of the virial ration T/U. However, we believe 
this effect to be negligible in a statistical study like ours. 

Our study shows that ~ 20 per cent of groups gener- 
ated by the FOF algorithm are not gravitationally bound 
when the Qa = 0.73 model i s ado pted. This result is in 
agreement with lDiaferio et al] (|l994h . who derived a similar 
result for the compact groups of galaxies. If we vary the ap- 
parent magnitude limit of the search from the original 13.2 
to 20.0, even more groups (~ 37 per cent) are unbound. 
This is not a negligible fraction considering that one widely 
accepted and applied method of calculating a group mass, 
from observations, is based on the assumption that groups 
found by the FOF algorithm are, in general, gravitationally 
bound systems. 

If we vary the value of the cosmological constant from 
the original 0.73 to 0.90, a slightly larger fraction of groups 
seems to be unbound when the apparent magnitude limit of 
13.2 is adopted. This result is intuitively reasonable. If the 
negative vacuum pressure of space is larger, gravitational 
force becomes 'weaker' and a smaller number of dark mat- 
ter haloes are formed and fewer groups are gravitationally 
bound objects. How does the fraction of gravitationally un- 
bound groups change when the negative vacuum pressure 
of space is lowered? If the value of the cosmological con- 
stant is put to 0.0, about the same fraction of groups (with 
Tnum = 13.2) are spurious as in the Qa = 0.90 cosmology. 
When the apparent magnitude limit is changed to 20.0, ~ 37 
per cent of the groups are spurious (for details, see Table (B]). 

In the low-resolution Qa = 0.73 simulation the frac- 
tion of gravitationally bound groups rises from 81.1 to 81.7 
per cent, when more recent values of the Schechter lumi- 



Table 6. Fractions of gravitationally bound groups of dark mat- 
ter haloes when different cosmological models and apparent mag- 
nitude limits have been adopted. 



Qa 




^ groups 


fbound {%) 


^isolated (%) 


0.00 


13.2 


2675 


79.2 


20.9 


0.00 


20.0 


6213 


64.5 


41.8 


0.73^ 


13.2 


1570 


70.1 


32.8 


0.73^ 


13.2 


1168 


81.1 


42.2 


0.73^ 


20.0 


6807 


62.7 


41.2 


0.90 


13.2 


238 


77.3 


53.1 


0.90 


20.0 


3678 


62.0 


36.5 



Note: Qa specifies the value of the cosmological constant ( = 
high-resolution and ^ = low-resolution simulation), rrin^ is the 
apparent magnitude limit of the search, Ngroups is the number of 
groups found from ten observation points, fbound is the fraction 
of gravitationally bound groups, and Ni^olated is the percentage 
of the isolated haloes which do not belong to any group or binary 
system. 

nosity function are adopted. Meanwhile the total number 
of groups decreases ~ 10.0 per cent. The fraction of grav- 
itationally bound groups rises from 81.1 to 82.8 per cent, 
when the values of the free parameters of Do = 0.37 Mpc 
and Vb = 200 km s"'^ are adopted. This result agrees with 
Frederic (1995a,b) who obtained similar result while study- 
ing group accuracy as a function of Do and Vb- Frederic 
(1995a,b) showed that smaller values of Do and Vb produces 
groups with greater accuracy and these groups should be 
gravitationally bound. 

Adopting the values of Do = 0.37 Mpc and Vo = 
200 km s"'^ have a significant effect to the total number 
of groups found from the simulations. The low-resolution 
Q.A = 0.73 simulation produces, in all, 1168 groups of dark 
matter haloes when the original values (Do = 0.63 Mpc 
and Vb = 400 km s~^) of the free parameters are adopted. 
When Do = 0.37 Mpc and Vb = 200 km s"^ are adopted, 
the total number of groups found from the low-resolution 
simulation drops to 661 while the fraction of isolated haloes 
rises from ~ 42.2 to ~ 64.9 per cent. Also, the group abun- 
dances changes significantly as the richest group found from 
the low-resolution simulation with Do = 0.37 Mpc and 
Vb ~ 200 km s~^ has only 15 members. These results are 
due to the fact that the limiting density enhancement of the 
search is inversely proportional to Dq. 

Our study shows that the 57a = 0.0 model produces 
about the same fraction of bound groups as the Q,a ~ 0.90 
model when the apparent magnitude limit is 13.2. However, 
the JIa = 0.73 model produces more gravitationally bound 
systems than the two other models we study when the orig- 



inal apparent magnitude limit is adopted. If the apparent 
magnitude limit is changed to 20.0, the Qa ~ 0.0 model 
produces a slightly larger fraction of gravitationally bound 
groups than the Qa = 0.73 or the Q.a ~ 0.90 simulations 
(for details see Table [6]). The use of the apparent magnitude 
limit of 20.0 means simply that every single dark matter 
halo in a simulation box is visible at the observation point. 
This result is not without bias as the simulation box is of fi- 
nite size and the edge effects might become significant, even 
using the periodic boundary conditions in the simulations. 

The calculation that determines whether a group is 
bound is based on three parameters: the total mass of the 
group, the relative velocity of the group members and the 
physical size of the group. In the following we will study 
how sensitive the result is on the values of these parame- 
ters. But first we will study the virial ratio as a function of 
the number of members in the group. Figure [6] shows the 
virial ratio T/U as a function of the number of haloes in 
the group, namely richness. The data come from the low- 
resolution ~ 0.73 simulation with mnm = 13.2. There 
are, in all, 1168 groups of haloes seen from ten different ob- 
servation points. The fractions of bound 'poor' and 'rich' 
groups are shown in Table [T] 

From Figure [S] we see that groups with more than 
10 members are most likely gravitationally bound and 
groups with three to five members are quite often unbound. 
iRameUa et all (|l99i1 ) argues that among groups with three 
members, 50 to 75 per cent of groups are spurious. They also 
conclude that for groups with more than three members the 
fraction of spurious groups is less than 30 per cent and may 
be as small as 10 per cent. Our findings are similar, and the 
fraction of bound groups w i th fou r or more me mbers is about 
as hig h as iRamella et al.l (119971 ) suggested. iRamella et al.l 
l|2002f) find that for groups with five or more members at 
least 80 per cent of the groups are probably physical sys- 
tems, but that 40 to 80 per cent of the groups with five or 
more members are bound groups. Our findings confirm the 
latter result. However, these results cannot be directly com- 
pared with ours, as slightly different values of the free pa- 
rameters are adopted for the FOF algorithm. Even-though 
the free parameters of the FOF seem to have only a small 
effect to the fraction of spurious groups in our study. 

Our findings for th e Qa = 0.73 cosmology are similar 
to lRameUa et all (1997') for 'poor' (three or four members) 
groups, as can be seen in Figure|S]and Table [T] However, our 
findings do not confirm the claim by IRamella et al.l (Il997l ) 
that among groups with three members, 50 to 75 per cent 
of groups are spurious as we find only ~ 23 per cent of 
groups with three members to be gravitationally unbound 
with miim ~ 13.2. For the apparent magnitude-limited sam- 
ple {miim = 13.2 comparable to HG82), we found ~ 21 per 
cent of the groups with three or four haloes to be spurious. 
'Rich' groups with five or more members are more often 
gravitationally bound than 'poor' groups, but the difference 
is relatively small at the apparent magnitude limit of 13.2 
as for 'rich' groups, we found ~ 15 per cent of the groups 
to be spurious. This is close to the upper limit proposed 
by Ramella et al. (1997, 2002). More details of our findings 
with different abundances and apparent magnitude limits 
are fisted in Table [T] 

In Figure [T] the virial ratio is plotted as a function of 
the velocity dispersion (t„ of the group. The plot shows a 
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Table 7. Fractions of gravitationally bound 'poor' and 'rich' 
groups of dark matter haloes generated from the low-resolution 
Qa = 0.73 simulation. 





^ haloes 


Ngroups 


fbound (%) 


13.2 


3 


448 


77.0 


13.2 


e [3,4] 


697 


78.5 


13.2 


> 4 


471 


84.9 


20.0 


3 


2786 


55.6 


20.0 


e [3,4] 


4147 


56.3 


20.0 


> 4 


2660 


72.6 



Note: rriHm is the apparent magnitude limit of a sample, N/^aloea 
is the number of haloes in a group, Ngroups is the number of 
groups found from ten observation points with appropriate num- 
ber of haloes, and fbound is the fraction of gravitationally bound 
groups. 
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Figure 6. Virial ratio (TU~^) versus the number of members in 
a group (Richness). Groups with more than 10 members are more 
often bound than 'poor' groups with three to five members. The 
data are from the low-resolution Q,a = 0.73 simulation when the 
apparent magnitude limit of 13.2 has been adopted. 

weak correlation in the sense that groups with large veloc- 
ity dispersion are more often gravitationally unbound than 
groups with small velocity dispersion. The linear correlation 
coefficient of 0.17 suggests that the correlation in FigureQis 
weak. However, the correlation is significant at level of 0.001 
and P{r) ~ 10'^ The rms line plotted in Figures [7] — [11] is 
of the form oc o-Jj or oc N{haloes)^ . The value of the 
parameter h of the rms line in Figure [7] is 6 = 0.10 ± 0.04. 

The weak trends are clearer if we scale the abscissa in 
both Figures [6] and [7| with the total mass of the group. Note 
that we use here the 'true' mass of a group rather than 
the 'observable' mass. Results are shown in Figures [8] and 
[9] More significant trends are now seen in both Figures, 
even-though the data are still scattered. The linear correla- 
tion coefficient of 0.32 suggests that a significant correlation 
exists between the number of haloes and the virial ratio 
when the first is scaled with the total mass of the group. 
The correlation in Figure [8] is significant at level of 0.001 
and P{r) < 10~^^. The slope of the rms line in Figure [8] is 
b = 0.82 ±0.03. 

Figure [9] shows a strong trend even-though the data are 
still quite scattered. The linear correlation coefficient of 0.40 
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Figure 7. Virial ratio {TU~^) versus the velocity dispersion (cr^) 
of a group (in km s~^). The rms straight hne has been fitted to the 
data. The data are from the low-resolution Qa =0.73 simulation 
when the apparent magnitude limit of 13.2 has been adopted. 




Figure 8. Virial ratio (TC/^^) versus the number of members 
(Richness) in a group when the latter has been scaled with the to- 
tal mass of the group (M^^o„p in A^q). The rms straight line 
has been fitted to the data. The data are from the low-resolution 
Qa = 0.73 simulation when the apparent magnitude limit of 13.2 
has been adopted. 
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Figure 9. Virial ratio {TU~^) versus the velocity dispersion ((t„) 
of a group (in km s~^) when the latter has been scaled with the 
total mass of the group (Mgroup in MQ).The rms straight line 
has been fitted to the data. The data are from the low-resolution 
Qa = 0.73 simulation, when the apparent magnitude limit of 13.2 
has been adopted. 




Figure 10. Virial ratio (TC/"^) versus the number of members 
in a group (Richness) when the latter has been scaled with the 
total mass of a group (MgJ'oup in A^q), and Qa = 0.0 and 
fniim = 20.0 has been adopted. The rms straight line has been 
fitted to the data. 



suggests that the correlation between the velocity dispersion 
of a group and the virial ratio of a group is quite strong. 
The correlation in Figure [9] is significant at level of 0.001 
and P{r) < 10~^^. The slope with standard errors of the 
rms line is now b — 0.90 ± 0.03. 

When the apparent magnitude limit is changed to 20.0, 
the trends of Figures [S] and |5] become stronger and the 
asymptotic standard errors for the rms lines become much 
smaller. The linear correlation coefficient of 0.62 shows that 
the correlation between the velocity dispersion of a group, 
and the virial ratio, is strong when the apparent magnitude 
limit of 20.0 is adopted. This correlation is significant at 
level of 0.001 and P(r) < 10"^^ What might be surprising 
is that changing the cosmological model, i.e. the value of 
the cosmological constant Qa, does not have a substantial 
influence on Figures 9 and 10. The number of groups found 
from different simulations varies a lot as a function of Q,a 
but the fraction of gravitationally bound groups do not (see 



Table |6j. For comparison with Figures [8] and [9] we show re- 
sults from the Qa ~ 0.0 simulation in Figures [TOI and [TT] In 
these figures the apparent magnitude limit of 20.0 has been 
adopted. Otherwise, Figures [10] and [TT] are comparable to 
Figures [8] and [9] The linear correlation coefficient in Figures 
llOl and fTTl is: 0.32 and 0.65, respectively. Both of the correla- 
tions are significant at a level of 0.001 and P{r) < 10~^^ for 
both samples. The asymptotic standard errors for the rms 
lines in Figures \M and [TT] are smaU: b = 0.83 ± 0.02 and 
b = 0.91 ± 0.01, respectively. 



7 DISCUSSION 

7.1 Probability functions of unbound groups 

In this subsection we briefiy discuss a method, which gives a 
theoretical probability of a group being gravitationally un- 
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Figure 11. Virial ratio (TfJ-^) versus the velocity dispersion 
(cr^) of a group (in km s~^) when the latter has been scaled with 
the total mass of a group (Mgroup in ^q), and = 0.0 
and mii^ = 20.0 has been adopted. Straight line is a rms fit to 
the data. 



Figure 12. Unbound probabiUty {P{TU-^ > 1.0)) versus veloc- 
ity dispersion ((Tu) of a group when the Q\ = 0.73 model, the 
apparent magnitude limit of 20.0, and a bin length of 0.2 are 
adopted. cr^M'^roup is in units of km s~^ (10^^h~^ Mq)~^ . 



bound. The mass of the groups is assumed to be known. 
In observations, estimations of group masses are less than 
precise at best, therefore the applicability of this method to 
observational data is merely hypothetical. The observable 
quantities we study are the velocity dispersion divided by 
the group mass o^Mg^oup and the mean pairwise separation 
divided by the group mass RpM~^„^^. 

To calculate the probability functions for the groups, 
the first step is to choose an appropriate bin length (gener- 
ally between 0.15 and 0.30) in the logarithm of the observ- 
able quantity. Then one calculates the number of groups and 
the number of gravitationally bound groups in each bin and 
divides the number of unbound groups with the total num- 
ber of groups in the bin. The logarithmic scale is chosen in 
order to lower the dispersion of the data and to assure a 
large enough number of groups in every bin. 

The probability functions for the velocity dispersion (t„ 
and the mean pairwise separation Rp, normalised to the 
group mass, are shown in Figures [T^ and lT^ Figure [T^] shows 
that at values larger than 1.8 (the horizontal dotted line in 
Figure \12\ it is more probable that the groups are gravi- 
tationally unbound when the Q.a = 0.73 model is adopted. 
The same result can also be inferred from Figured but with 
lower confidence. 

The f^A = 0.0 simulation gives a probability func- 
tion which is comparable to the probability function of the 
Q.A = 0.73 model. It shows a similar linear growth as the 
probability function of the = 0.73 model. However, the 
function is shifted along the horizontal axis. This shifting 
originates from the variation of the group masses and ve- 
locity dispersions (see Figure [Sjl. The change of the appar- 
ent magnitude limit does not have any significant effect on 
Figure 1121 When the apparent magnitude limit of 13.2 is 
adopted, a smaller number of haloes and groups are ob- 
served, which enlarges the variations between bins and gives 
a worse fit to a straight line. 

The quantity a"^Mg^^o„p, studied in Figure 1121 has a 
loose connection to the kinetic energy T. This connection 
explains the fact that the lower normalized values of the 
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Figure 13. Unbound probability (P{TU~^ > 1.0)) versus mean 
pairwise separation (Rp) when the £7a = 0.73 model, the apparent 
magnitude limit of 20.0, and a bin length of 0.25 are adopted. 
-RpMgroup is in units of Mpc {lO^^h''^ Mq)-'^. 



a higher probability, and larger normalized values, loosely 
meaning the larger kinetic energies, gives unbound groups 
with a higher probability. 

The probability function of the mean pairwise separa- 
tion (Figure I13|) shows a similar linear growth as the prob- 
ability function of the velocity dispersion in Figure [121 The 
variation from bin to bin is somewhat larger in the probabil- 
ity function of the mean pairwise separation due to the fact 
that the mean pairwise separation is not strictly the size of 
the group but it includes projection effects. The quantity 

Rp ^^group 5 



quantity cr„Mg 



gives gravitationally bound groups with 



^ studied in Figure 1131 is inversely proportional 
to the potential energy U if the mean pairwise separation 
is identified as the real size of a group. The difference be- 
tween Figures \n\ and [THl is understandable, as the velocity 
dispersion and the mean pairwise separation are not strictly 
connected to each other, although some loose relation exists 
as the equation of the mean pairwise separation includes the 
mean group radial velocity. 

The Oa = 0.90 simulation show similar probability 
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functions as the = 0.0 and the SI a =0.73 simulations. 
Figures [12] and [13] show that the f^A does not have any sig- 
nificant effect for the fraction of gravitationally unbound 
groups. This result can also be inferred from Table [6] The 
small ef fect is hardly surp rising as the theoretical studies 
(see e.g. iLahav et al.lll99"ll ) have predicted that the Q,a has 
little effect on the dynamics at the present epoch. 



8 CONCLUSIONS 

We have shown that the ACDM cosmology can produce 
groups of dark matter haloes comparable to observations of 
groups of galaxie s when the FOF algorithm based on that of 
iHuchra fc Gelled (|l982l ) is adopted and the dynamical prop- 
erties of groups are studied. Our groups from cosmological 
simulations are, in general, in a moderate agreement with 
observations, although a straight K-S test fails in most cases. 
Our Qa = 0.73 simulations are in satisfactory agreement 
with observations as the number densities of group proper- 
ties are usually within 2a errors, or less, from the HG82 and 
the UZC-SSRS2 group abundances. The agreement between 
simulations and observations is good when the velocity dis- 
persion and the 'observable' mass of groups is considered. In 
these cases the applied K-S test is approved when the high 
resolution Qa = 0.73 simulation is considered. The moder- 
ate agreement between simulations and observational data 
suggests that gravitational force alone is sufficient in order 
to explain the dynamical properties of groups of galaxies. 

We have also shown that in general about 20 per cent 
of the groups of haloes generated with the algorithm pre- 
sented in the HG82 are not gravitationally bound objects. 
The fraction of gravitationally bound groups of dark matter 
haloes varies with different values of the apparent magni- 
tude limits. When the apparent magnitude limit is raised 
from the original 13.2 to 20.0, a larger number of spurious 
groups are found. The larger fraction of unbound groups 
with miim = 20.0 could be explained by the fact that 
more interlopers are included into groups, when the appar- 
ent magnitude limit is increased. However, this analysis is 
beyond the scope of this study. 

In general, a larger number of 'rich' groups are found 
when the apparent magnitude limit is lowered. This origi- 
nates from the fact that more light haloes at close proximity 
to more massive haloes become visible and those light haloes 
are included into the groups. When the magnitude limit is 
raised from the original value of 13.2 to 12.0, a slightly larger 
fraction of the groups are found to be gravitationally bound. 
In general, fewer groups (in absolutely number) are found 
and these groups are 'poorer'. 

Small differences are found when the fractions of gravi- 
tationally bound 'poor' and 'rich' groups are studied. 'Rich' 
groups with more than four members are more often gravita- 
tionally bound tha n 'poorer' groups. Th is result agrees with 
previous ones (e.g. iRamella et al.l 20021). O ur results do not 
confirm the claim bv iRamella et al' ( ^93) who argued that 



50 to 75 per cent of groups with three members are spurious. 
Our results show that ~ 77 per cent of groups containing 
only three members are gravitationally bound when the ap- 
parent magnitude limit of 13.2 is adopted. 

When the value of the cosmological constant JIa is var- 
ied, the fractions of unbound groups change only slightly. 



This is somewhat surprising as it would be intuitively ex- 
pected that a larger value of the dark energy would lead 
to a greater number of groups that are not gravitationally 
bound. Some variation is observed when the fraction of grav- 
itationally bound groups is studied as a function of the cos- 
mological constant, but, in general, a significant number of 
groups remains unbound in all cases of Qa- 

When the values of the free parameters of the FOF 
algorithm are varied, the fraction of gravitationally bound 
groups can be raised from ~ 81 to ~ 83 per cent. A greater 
difference is observed when the fraction of isolated haloes 
is studied. Varying the values of Do and Vo makes a great 
difference, raising the fraction of isolated haloes from ~ 42 
to ~ 65 per cent. In general, we do not find any significant 
difference in the fractions of gravitationally bound groups 
when different values of Do and Vo or parameters of the 
Schechter luminosity function are adopted. 

In observations, the crossing time of a group is often 
taken as an indicator of the virialization. We do not find any 
correlation between the virial ratio and the crossing time of 
a group. This result does not depend on the chosen value of 
the apparent magnitude limit of the search, or the cosmo- 
logical model adopted. The lack of the correlation between 
these two variables calls into question the crossing time as 
an estimator of the virialization. 
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